A traveling wave decelerator for neutral polar molecules 
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Recently, a decelerator for neutral polar molecules has been presented that operates on the basis 
of macroscopic, three-dimensional, traveling electrostatic traps (Osterwalder et al., Phys. Rev. A 
81, 051401 (2010)). In the present paper, a complete description of this decelerator is given, with 
emphasis on the electronics and the mechanical design. Experimental results showing the transverse 
velocity distributions of guided molecules are shown and compared to trajectory simulations. An 
assessment of non-adiabatic losses is made by comparing the deceleration signals from ^'^CO with 
those from ^^CO and with simulated signals. 
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I. INTRODUCTION 

Interest in the production of cold molecules stems from 
a variety of research areas, most notably collision dynam- 
ics, the study of long-range interactions, quantunr com- 
putation, and high-resolution spectroscopy (for a recent 
overview see the individual chapters of reference 1). 

Two principal approaches can be distinguished in the 
production of cold molecules: either the molecules are as- 
sembled from ultracold atoms using electric fields, mag- 
netic fields, or lasers, or one starts from warm or fast 
sanrples and brings them to low average translational 
temperatures by controlling the velocityjii^ Variants of 
the latter approach include using collisions with a cold 
buffer gas? or between molecules in crossed-beams,"' by 
mechanical means such as rotating nozzles^i^ or through 
the controlled interaction with electric^Ti^ magnetic;^""— 
or optical fields. ^^ 

One sinrple way to use electric fields to passively con- 
trol a nrolecular velocity distribution is by filtering out 
fast molecules using the competition between an electro- 
static force and a centrifugal force in a bent guideJ^"— 
This does not change the velocity of the individual 
molecules but merely extracts the slow ones from a ther- 
mal distribution. Active control of the forward velocity 
of molecules with electric fields was first demonstrated 
in 1999}i^ rapidly followed by the first trapping of neu- 
tral molecules jiiii^ Further developments of this technol- 
ogy have enabled the deceleration of large molecules in 
high-field-seeking stateSjiS, the trapping of molecules in 
static and ac-fields /^i^° and the storage of molecules in 
a synchrotron.^^'-- 

For most of the studies mentioned above, consider- 
able density both in velocity space and in position space 
are of great inrportance. A significant developinent in 
Stark deceleration was thus the recognition that such de- 
vices should be operated in a higher-order mode in order 
to suppress the coupling of longitudinal and transverse 
oscillations and the losses following from them i^^i^^ Re- 
cently, a new approach has been presented to completely 
avoid this problem, namely the use of traveling three- 



dimensional traps to capture and decelerate moleculesj^ 
This approach builds on ideas implenrented in a miniatur- 
ized Stark decelerator on a chip)^"— and a similar prin- 
ciple has also been applied to the deceleration of Rydberg 
states^Si. In contrast to a conventional Stark decelerator, 
where effective moving traps are obtained by periodic 
switching of the electric fields between two static config- 
urations, these devices employ continuous modulation of 
the voltages to produce real, continuously nroving three- 
dimensional traps. In the chip decelerator, this modula- 
tion is chosen to produce a chain of traps above a periodic 
structure of electrodes, which in this case are stripes on 
a dielectric substrate. The modulation frequency deter- 
mines, for a particular electrode geometry, the velocity 
of the traps and can be chirped to produce acceleration 
or deceleration. In the present decelerator, the stripes of 
the chip decelerator are replaced with three-dimensional 
ring electrodes.^^ The decelerator is scaled up by a fac- 
tor of 100 in comparison with the chip decelerator, which 
results in dimensions comparable to conventional Stark 
decelerators. Based on trajectory calculations, the six- 
dimensional phase space acceptance for guiding of OH 
molecules at a constant velocity is determined to be a fac- 
tor of three larger than in the state-of-the-art switched 
Stark decelerator f24 and during deceleration, the same 
relative phase space acceptance is maintained for an ac- 
celeration which is a factor of three larger. Furthermore, 
the molecules are loaded into the traps at high velocity, 
and the loading of stationary traps is simply bringing the 
traveling traps to standstill. 

The present paper gives a detailed account of the de- 
sign and construction of this new decelerator and shows 
experimental data together with the corresponding tra- 
jectory simulations to explain the operation principle of 
the device and the dynamics of individual molecules in- 
side the decelerator. The paper is structured as follows: 
section |ll] describes the general principle of operation 
pi Ap . the mechanical design (jIIB[) . and the required 
electronics (jIICp . Section IIIBI includes the particular 
approach chosen to assemble the structure. Section IIIII 
presents experimental results that characterize the appa- 



ratus, and theoretical calculations to support the inter- 
pretation of these measurements. The paper is summa- 
rized in section ITVl where also an outlook is given, along 
with future applications of this design. 



II. EXPERIMENTAL 

The following sections give a general overview on the 
operation principle of the decelerator and describe in de- 
tail how the device is designed and constructed and how 
the required electronics are built. 



A. General considerations 

Applying a spatial sine-modulated electric potential to 
a periodic array of ring electrodes produces an electric 
field distribution on the inside of the array with two on- 
axis minima in the field magnitude per period (see figure 
[T]). Due to the cylindrical symmetry, these minima rep- 
resent three-dimensional traps for polar molecules. 

The required potential on the individual electrodes can 
be expressed as 



Vn{t) = Va sin(-0o(t) + 2Trn/N) 



(1) 



with the waveform amplitude Vq and the time-dependent 
phase offset (f>o{t). The electrode-dependent shift is de- 
termined through the periodicity N of the specific array. 
A linear increase of (j)o{t) with time continuously shifts 
the traps along the cylinder axis. The time-dependence 
of the phase offset can be expressed as 



^W 



2tt I v{T)dT, 
/o 



(2) 



where the frequency v{t) is the modulation frequency 
of the voltage on each individual electrode. Since one 
oscillation of the waveform moves the trap over one pe- 
riod, the velocity of the traps is given by Lv{t), where 
L is the period length. The time-dependence of the fre- 
quency leads to acceleration or deceleration of the traps. 
After calculating the acceptance for many different geo- 
metrical parameters, a suitable combination of these was 
found with 4 mm inner electrode diameter, 0.6 mm elec- 
trode thickness, 12 mm period length, and periodicity 
8. Lower periodicity means fewer points to sample the 
spatial sine, and thus less smooth electric field distribu- 
tions. Higher periodicity improves the electric fields but 
makes the setup technically more challenging. The field 
distribution produced by a larger number of electrodes 
per period alone is not necessarily better: in optimizing 
the number of electrodes, it was found that - while keep- 
ing all other parameters the same - a decelerator with 8 
electrodes per period performs better than one with 10, 
though a decelerator with 6 is significantly worse than 
either of these. 
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FIG. 1: Calculated electric fields for the trap located in the 
plane of an electrode (A and B), and between two electrodes 
(C and D). Panels A and C show the radial dependence of 
the electric field magnitude, panels B and D show contour 
plots of the electric field magnitude. Line types in panels A 
and C correspond to cross sections at the positions indicated 
by lines of the same type in panels B and D. Contours start 
at 5 kV/cm around the center of the trap and are spaced 
by 5 kV/cm. The thick lines represent 10 kV/cm and 20 
kV/cm. Gray circles on the right of panels B and D show 
the positions of the electrodes. Numbers next to the circles 
indicate the voltages applied to the individual electrodes for 
both cases. Panel E shows the longitudinal dependence of 
the electric field magnitude for the cases from panels B (solid 
line) and D (dashed line). 



In order to create sufficiently deep traps moving at the 
required velocity for this geometry, the electronics must 
produce waveforms with Vo wlO kV, with a tunable fre- 
quency up to sa25 kHz (corresponding to a velocity of 300 
m/s). The phase of the waveforms between neighboring 
electrodes is shifted by 7r/4. 

Contour plots for single traps resulting from Vb=10 kV 
sine potentials used with this design are shown in figures 
[T)3 and D for the cases where the trap is positioned in the 
plane of an electrode (corresponding to 0o = 0) and be- 
tween two electrodes (0o = ""/S), respectively. Contours 



are spaced by 5 kV/cm, starting at 5 kV/cm. Thick lines 
show the contours for 10 kV/cm and 20 kV/cm, respec- 
tively. The gray circles to the right of panels B and D 
show the positions of the electrodes. Numbers next to 
the circles show the voltages that are applied to each 
electrode in the two cases. Panels A and C show the 
radial dependence of the electric field magnitude at po- 
sitions designated by the lines with corresponding line 
types in panels B and D. Each line represents the cross 
section for different distances from the trap minimum. In 
both cases, the solid, dotted, dash-dotted, long-dashed, 
and long-short-dashed lines represent mm, 0.75 mm, 
1.5 mm, 2.25 mm, and 3 mm from the trap minimum; 3 
mm from the trap minimum is the position of the sad- 
dle point toward the next trap. The on-axis longitudinal 
potential is shown in panel E for the two cases. Here, 
the solid (dashed) line represents the case with (/)o = 
{(j)a — 7r/8). Panel E shows that, longitudinally, the po- 
tential does not change when moving the trap from in- 
plane to between two electrode planes. This is also con- 
firmed through the height of the saddle point in panels 
A and C. Longitudinally, the traps are almost 40 kV/cm 
deep when sine waves with amplitudes of 10 kV are ap- 
plied to the electrodes. For a linearly Stark-shifted state, 
the trap has the shape of | sin(a;)|, which is nearly linear 
close to the minimum. While the shape of the longi- 
tudinal potential is very stable and independent of (Jjq, 
the radial potential does undergo some oscillations. The 
trap is almost 50 kV/cm deep when located in the plane 
of an electrode, but sal3 kV/cm when between two elec- 
trodes. For the a^lli{v — 0, J — 1) state of CO used 
here this leads to a transverse acceptance of ss ±22 m/s 
and ~ ±11 m/s, respectively. Trajectory simulations on 
a virtual decelerator with 10 m length demonstrated that 
the effective transverse acceptance lies in between these 
values and is close to the one obtained when averaging 
the radial potential along the guide axis. Longitudinally, 
the acceptance for guiding at constant velocity is « ±19 
m/s^ 



B. Mechanical design and construction 

The upper part of figure [5] is a photograph of the 
top end of the decelerator. The central part shows the 
electrodes, visible as horizontal bars, and the electrode 
mounting rods, as described below. The lower part of the 
photograph shows the mounting structure with the main 
octagonal stainless-steel tube and bronze plates that hold 
the electrode mounts. An overview of the decelerator de- 
sign, including its mount, is shown in figure [2]A_. Figures 
[33-D show, respectively, a top view of the decelerator, 
a single electrode mount, and an electrode with a cross 
section through the mount. The decelerator is positioned 
vertically inside a high- vacuum chamber. Molecules are 
produced in a pulsed supersonic expansion and enter the 
decelerator chamber through a skimmer, visible at the 
bottom of figure ^^, that is placed as close as possi- 




FIG. 2: Upper section: Photograph showing the upper end of 
the decelerator structure. Bottom section: Mechanical design 
of the traveling wave decelerator. A) Cross section through 
the entire decelerator. The arrow labeled "CCD" shows the 
position of the CCD camera used for the phosphorescence 
imaging. B) Top view of the decelerator and mount structure. 
C) Single rod containing one set of 40 ring electrodes. D) 
Top view of a single electrode with cross section through the 
mounting rod. The active element of the electrode is the 
circular section at the top in this figure. 



ble to the beginning of the decelerator. In the case of 
the CO experiments presented below, a UV laser crosses 
the expansion directly before the skimmer to excite the 
CO molecules to the upper A-doublet component of the 
metastable a^ni(ti = 0, J = 1) state. This particular 
state has a phosphorescence lifetime of 2.63 ms.'^'' 

The decelerator is formed from 320 ring electrodes that 
are mounted on 8 separate 6 mm-diameter stainless steel 
rods. It is worth remembering that the number of elec- 
trodes here can not be compared to the number of decel- 
eration stages in a traditional Stark decelerator; the fig- 



ure of merit is instead the electric field gradient that can 
be generated, and this determines the length of the de- 
celerator, which in turn defines the number of electrodes. 
The 8 rods are each mounted on ceramic posts that are 
attached via adjustable bronze bars to the main octag- 
onal mounting structure (see figure [2j3). The mechani- 
cal point of reference of the entire setup is the center of 
the base flange that holds the skimmer and supports the 
decelerator. This construction, and the vertical orienta- 
tion of the decelerator, greatly simplifies the alignment 
by making it independent of the main vacuum chamber. 
A magnification of one of the rods is shown in panel [5p. 
The individual electrodes are spaced by f2 mm which 
defines the period length. All 8 rods, hand-picked for 
straightness, are identical and are mounted with a lon- 
gitudinal offset of 1.5 mm from their nearest neighbors, 
a value which is given by the required center-to-center 
distance of neighboring electrodes. The arrangement of 
successive electrode rods is counter-clockwise when seen 
from the skimmer, forming a left-handed helical pattern. 
A single electrode is shown in panel [2p. The active 
part of these electrodes is the top, circular part with an 
inner diameter of 4 mm. The bottom part is the handle 
via which the electrode is mounted to the main rod. Each 
electrode is formed from 0.6 mm diameter Ta wire in a 
three-step procedure: 

1. The wire is stretched, polished, and cut to the ap- 
proximate length required for the electrode; 

2. Each piece is bent into a rectangular "U"-shape by 
forcing it around the corners of a 13 mm (« 7r-4 
mm, found by trial and error) hardened steel plate; 

3. In order to obtain the inner circular shape, each 
electrode is bent around a 4 mm stainless steel rod. 

Tantalum wire is used because of its high resilience to 
electrical discharges, and because it can be bent into the 
desired shape quite easily. '^^ The stretching of the wire is 
necessary to ensure straight wires, but it also reduces the 
diameter to 0.57 mm. Step 2 is crucial because it ensures 
that the wires are bent at right angles at the point where 
the two ends of the circular bend meet. All attempts 
to bend the electrodes directly into this shape, without 
the pre-bending, resulted in electrodes with cross sec- 
tions that resembled teardrops rather than circles. Con- 
siderable deviation from the overall cylindrical symme- 
try of the decelerator is obtained if the individual elec- 
trodes are not as close to circular as possible. The above 
procedure guarantees that all 320 electrodes are identi- 
cal to within high accuracy, and that the electrodes are 
produced without being scratched, thus rendering them 
more resistant to discharges when high voltages are ap- 
plied. The electrodes are then cleaned in an ultrasonic 
bath and mounted on the main electrode rods. Clean- 
ing is an important step in order to remove traces of the 
polishing material from the surface of the electrodes, and 
this step significantly enhances the resilience of the elec- 
trodes towards electrical discharges. 



Each of the mounting rods contains 40 holes for the 
electrodes, one of which is seen in figure [5p, and two 
tapped holes for the mounting. The ring electrodes are 
pushed into the holes and pre-aligned using a PVC block 
that fixes the distance between the rod surface and the 
electrode bend. The Ta electrodes are then fixed to the 
steel rods by pushing one or two pieces of bent, 0.3 mm 
piano wire into the spaces left between the Ta wire and 
the walls of the hole. This ensures that the electrodes do 
not slip out, but leaves the required flexibility for final 
alignment. At this stage, the electrodes are positioned to 
within saO.5 mm, insufficient for the present purposes. 

The exact positioning of the electrodes on the indi- 
vidual rods is complicated by the fact that their posi- 
tions have to be measured without mechanical contact, 
both in order to avoid damage on the surface of the elec- 
trodes, and because any contact might bend the wires. 
A solution is obtained by pre-mounting the rod with the 
ceramic mounts and the bronze bar (in the same combi- 
nation as in the final setup) on an aluminum block. This 
block is then positioned on a digital microscope that is 
built by fixing a CCD camera equipped with a macro- 
zoom lens to the head of a mill. The aluminum block, 
together with the pre-mounted electrode assembly, is ac- 
curately positioned on the moving table of the mill. In 
this way the relative positions of individual electrodes on 
each rod can be measured and set with a limit given by 
the resolution of the CCD. By placing the mount with the 
electrodes horizontal or vertical, it is possible to align and 
position the electrodes relative to the main rod and rel- 
ative to each other with an estimated accuracy of ±0.02 
mm. 

Finally, the 8 rods are mounted on the main cylinder 
and aligned to each other. This last step is done by po- 
sitioning one CCD camera above and one to the side of 
the decelerator, and a theodolite below it. Using custom- 
made alignment tools, each rod is positioned relative to 
its neighboring rods and to the main axis of the decel- 
erator. The vertical alignment can only be measured at 
one particular point, but the very accurate positioning 
of the individual electrodes on the mounting rod ensures 
the same degree of alignment throughout the entire struc- 
ture. 



Parameter 



Electrode diameter 

Distance between adjacent electrodes 

Electrode center to mounting rod 

Electrode offset from beam axis 

Electrode tilt around handle 

Rod in plane perpendicular to beam axis 

Rod position along beam axis 

Tilt of rod perpendicular to beam axis 



Precision 



±0.02 mm 
±0.02 mm 
±0.02 mm 
±0.02 mm 
±10 nirad 
±0.03 mm 
±0.02 mm 
±2 mrad 



TABLE I: List of estimated precision for critical parameters 
of the ring electrode decelerator. 



A quantitative assessment of the overall precision in 
the alignment is difficult. Table |T] summarizes the main 
estimated deviations from the ideal structure. These esti- 
mates base on a critical evaluation of each of the produc- 
tion steps, as well as on test measurements on selected 
electrodes and electrode positions after the alignment. 
The electrode diameter is defined through the bending 
around the 4 mm rod, and the deviation stated in the 
table has to be read as possible deviations from a per- 
fectly round shape, as well as a deviation from the ideal 
diameter. Due to the use of optical tools in the alignment 
procedure, the most difficult, and thus also least precise, 
part is the accurate positioning of the main electrode 
mounting rods in the plane perpendicular to the guide 
axis, mainly because of the mutual visual obstruction of 
the electrodes. 

An attempt to characterize the effects of these de- 
viations has been made by simulating the resulting 
electric fields in three dimensions using finite element 
calculations i^ A horizontal displacement of a single elec- 
trode moves the trap off-center when it is located in 
the plane of that electrode and has little effect other- 
wise. Displacement of an entire rod thus leads to a peri- 
odic jumping of the traps, with a frequency given by the 
modulation frequency. Vertical displacement changes the 
electric field between two electrodes but has little effect 
on the field inside of the decelerator. Longitudinal dis- 
placement of an entire rod can thus lead to a periodic de- 
formation of the traps to produce harder or softer walls 
with a periodicity also determined by the modulation fre- 
quency. While the modulation frequencies in the exper- 
iments done so far were substantially higher than the 
oscillation frequencies of the molecules inside the traps, 
which are in the range of 1 kHz, it still has to be as- 
sumed that the mechanical inaccuracies to some extent 
contribute to the increased losses observed experimen- 
tally when compared with the theoretical predictions. In 
the case of longitudinal displacement, a higher risk is 
however the increased danger of electrical discharges. 
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FIG. 3: Simplified circuit diagram for the generation of the 
high- voltage waveforms. AWG designates the arbitrary wave- 
form generator. Each of the eight channels in the current 
setup was driven by a separate amplification line fike this. 



considering the magnitude of other error sources, this 
effect can be neglected. For every channel, there is a sep- 
arate amplifier (Servowatt DCP 780/60 HSR) that raises 
the waveform amplitude from ±0.8 V to ±40 V. The high 
power of each amplifier (750 W) is required to drive a high 
current (±25 A) through the primary winding of a trans- 
former (custom made by Weiers und Partner) to produce 
the final ±10 kV waveform. The phase of the individ- 
ual channels is conserved during the amplification. A 
slight frequency dependence of the overall-amplification 
is, however, observed. In order to avoid excessive fluctu- 
ation of the trap shape, the amplitude on each electrode 
mounting rod is first checked as a function of frequency. 
The observed deviations from a constant amplitude are 
then corrected by modifying the input waveform. Since 
neighboring electrodes capacitively influence each other, 
this procedure has to be done iteratively to ensure con- 
stant and stable amplitudes. The process converges fast, 
however, and not more than two iterations were neces- 
sary in general to produce waveform amplitudes constant 
to within a few percent. 



C. High voltage- waveform generation 

A diagram of the electronics used to drive a single 
electrode rod is shown in figure [31 Computer-controlled 
waveforms with amplitudes of up to ±1 V are initially 
generated in an arbitrary waveform generator (AWG; 
Wuntronic, DA8150) for all eight channels simultane- 
ously. Programming of the AWG gives full control over 
the phase and amplitude of each of the waveforms at ev- 
ery moment during the deceleration, and for each of the 
channels separately. The AWG provides 150 MS/s at 12 
bit vertical resolution which corresponds to w6 kS per 
period at 25 kHz modulation. The 0.5 mV vertical reso- 
lution corresponds to only 5 V in the final waveform. The 
amplification procedure filters out high frequencies, and 
limited vertical resolution produces a slight uncertainty 
in the final voltages applied to the electrodes. However, 



III. RESULTS 

The operation of the apparatus is demonstrated in the 
following sections. In the first section, the phosphores- 
cence of the metastable CO molecules is used to image 
their motion along the decelerator axis and the forma- 
tion of density gradients near the trap positions. In the 
second section, the transverse distributions are investi- 
gated with application of dc voltages to the electrodes, 
thus forming an electrostatic guide. The third section de- 
scribes the results of the imaging of molecules that were 
guided in modulated electric fields. Finally the effect of 
potential non-adiabatic transitions in CO is studied by 
comparing arrival time distributions of ^^CO and ^^CO. 
These results complement the measurements presented 
previously. ■^^ There, arrival time distributions for -'^'^CO 
were presented for guiding with modulated potentials and 
for deceleration. 
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FIG. 4: Phosphorescence images of CO packets flying through 
the decelerator structure, recorded with a CCD camera from 
the side of the decelerator chamber. The central part of panel 
A shows an image recorded with a delay of 617 ps from the 
laser pulse and a width of 2 jis. The delay is chosen to cap- 
ture the trap exactly between two electrodes. The lower part 
of panel A shows the electric field magnitude at that time, 
and the curve at the top is a vertical integration of the image, 
showing the spatial dependence of the molecule density. The 
intensity scale is cut at «20% of the peak maximum to em- 
phasize the smaller peaks on the side (see text) . The zeroes in 
the curve are the shadows of the electrodes. Panel B shows a 
series of images, starting at a delay of 525.5 ps after the laser 
and with 5 lis wide gates. The time delay between subsequent 
images is 20 ps, corresponding to 1/2 period at 25 kHz. Above 
and below the images are plots of the electric field magnitude 
which has the same shape for each image. Vertical dashed 
lines show the calculated positions of the traps to guide the 
eye. 



Phosphorescence images of molecules in the 
decelerator 



The metastable state of CO used here has a lifetime of 
2.63 ms,^° and it decays to the electronic ground state 
by emission of a photon with wavelength between 206 nm 
and 240 nm<^ This phosphorescence can be imaged, us- 
ing an image-intensified CCD camera, to follow the den- 
sity of CO molecules as a function of time and position. 

Figure S] shows images taken from the side of the 
decelerator about 15 cm from its entrance (see figure 
[2]A) . These images were taken through the gaps between 
the decelerator electrodes, and the shadows of the elec- 
trodes are visible as a periodic intensity modulation. The 



molecules were guided through the decelerator at a con- 
stant velocity of 300 m/s. Figure I31A. shows a single 2 ys- 
long snapshot recorded at a fixed delay of 617 ps rela- 
tive to the laser pulse that produces the metastable CO 
molecules. This delay was chosen such that the densest 
part of the distribution would be located exactly between 
two electrodes. The horizontal axis corresponds to the 
molecular beam axis, with the source on the left and the 
detector on the right. The graph above the image shows 
the same data integrated over the transverse coordinate. 
The intense region at the center of the image is produced 
by the single trap that is filled by the molecular beam; 
the false color image of the electric field magnitude at 
the bottom of the panel shows the position of the traps 
(darker areas) at the time of the snap-shot. On either side 
of the main peak are molecules that jumped over the bar- 
rier to the previous or next trap. Because of their energy, 
they continue to pass over all saddle points along the de- 
celerator. Since they slow down on the saddle points, 
there is a slight increase in the density which is visible 6 
electrode spacings away on either side of the main peak. 
The structure becomes more pronounced by the time 
the molecules have passed through the entire decelerator, 
and they are very clear in the arrival time distributions 
shown, e.g., in figure [T] The detection efficiency of the 
imaging is relatively low. The 5 cm-diameter lens 35 cm 
from the center of the decelerator only collects 0.13% of 
the light emitted by the phosphorescing molecules, and 
the imaging system itself is only about 10% efficient. Ad- 
ditionally, during a 2 ps exposure, only about 0.07% of 
the molecules will phosphoresce. The resulting, overall 
detection efficiency is roughly 10~^. From this and the 
total number of photons observed in the images, we can 
infer that about 10^ molecules are guided in the main 
trap at a density of approximately 10^ molecules/cm"^. 



By changing the time at which the image is taken, the 
motion of the cloud of molecules through the decelera- 
tor can be followed. Panel B of figure H] shows a series 
of 5 ps-wide images taken at 20 ps intervals, starting at 
525.5 ps. A guiding velocity of 300 m/s requires 25 kHz 
waveforms, and the time interval corresponds exactly to 
one half period. Since there are two field minima per 
period, the electric field distribution (shown at the top 
and bottom of the panel) is the same for each image. The 
dashed vertical lines show the positions of the traps at the 
respective times. The images show a cloud of molecules 
that is confined both transversely and longitudinally as it 
is guided through the device in a single moving trap. The 
molecules outside this trap are transversely confined, but 
will spread out longitudinally as they proceed through 
the decelerator; the only remaining structure in the ar- 
rival time distribution is given by the bunches formed on 
the saddle points, as described above. 
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FIG. 5: Transverse position distribution of the molecules after 
the decelerator with kV, ±1 kV, ±2 kV, and ±3 kV applied 
to the electrodes. Panel A: electric field distribution in the 
decelerator with alternating positive and negative voltages on 
the electrodes. Panel B: two-dimensional transverse distribu- 
tions at the detector for the various applied potentials. Panel 
C: angle- integrated radial intensity profile of the distributions 
shown in panel B. 



B. Guiding of molecules with static fields 

Errors in the operation of the decelerator can stem 
from imperfections in the construction and in the wave- 
forms. By applying static potentials to the electrodes, 
the errors due to imperfections in the construction can 
be examined independent of the waveform errors. If al- 
ternating positive and negative potentials are applied to 
the electrodes, the resulting distribution of electric field 
magnitude has a minimum on axis (shown in panel A of 
figurein]) and increases steeply near the walls of the tube. 
This electric field distribution generates a guide that con- 
fines the molecules in low field seeking states transversely 
only. 

Although there is also a longitudinal modulation of the 
electric field strength (including a point of zero electric 
field at the center of each electrode), this modulation 
is too weak to significantly affect the longitudinal mo- 
tion of the molecules; application of guiding potentials 
increases the flux without significantly changing the ar- 
rival time distribution. The transverse position distribu- 
tion several centimeters behind the decelerator, on the 
other hand, shows significant structure that can be used 
to characterize the decelerator. To measure this distri- 



bution, an imaging detector (microchannel plate (MCP) 
detector, coupled to a phosphor screen) is placed 18 cm 
behind the last electrode of the decelerator. Metastable 
CO molecules that hit the MCP generate an Auger elec- 
tron which is subsequently amplified by the MCP stack. 
The resulting spot on the phosphor screen is recorded 
on a CCD camera. The center position of each spot is 
determined using a centroiding algorithm. By switch- 
ing the voltage applied to the front of the MCP stack, 
detection is restricted to molecules with an arrival time 
corresponding to a longitudinal velocity of 300±6 m/s. 

The images in figure[Sj3 show the position distributions 
of molecules at the MCP with kV, ±1 kV, ±2 kV, or 
±3 kV applied to the electrodes of the decelerator. At 
zero volts (upper left image in figure [5}3), the molecules 
form a circular distribution, the size of which is mainly 
determined by the geometric aperture formed by the di- 
ameter and length of the guide, and the distance to the 
imaging detector. Two main changes are observed as 
the potential on the electrodes increases: 1. The over- 
all intensity increases, and 2. in addition to the intense 
central disk, larger rings appear that grow with increas- 
ing voltage. As the potential is increased, molecules with 
larger transverse velocities can be guided through, result- 
ing in the increased overall intensity. These molecules 
also execute more oscillations in the device, producing 
more rings that grow in diameter. While the images are 
almost circularly symmetric, some asymmetries are also 
visible. The prominent left-right asymmetry, particularly 
in the outer rings, results in part from an off-center place- 
ment of the laser excitation at the beginning of the ex- 
periment, and it can be influenced by moving the laser. 
An up-down asymmetry of a similar degree is obtained 
by slightly detuning the laser. This detuning leads to the 
selective excitation of molecules with a velocity compo- 
nent along the laser beam that are Doppler-shifted into 
resonance. The bandwidth of the laser corresponds to a 
velocity range of w30 m/s which is in a range similar to 
the acceptance of the guide. A subtle eight-fold symme- 
try is also visible in some of the rings, which results from 
the eight different positions of the electrode indentations. 
These result from the shape that was chosen for the ring 
electrodes here (see figure[5p). The subtlety of this struc- 
ture in the images indicates that the indentations do not 
significantly affect the operation of the decelerator. 

The left (right) panel of figure [Sp shows the experi- 
mental (calculated) angle-integrated radial intensity pro- 
file of the images in panel B. The calculations were done 
via two-dimensional trajectory simulations that assume 
cylindrical symmetry. In all cases except with no volt- 
age applied, the relative number of molecules reaching 
the detector is accurately reproduced by the trajectory 
simulations. In the field free case, the lower number of 
molecules reaching the detector in the simulations can be 
qualitatively explained by the existence of non-low field 
seeking states; these states are present in the experiment 
and travel along the same straight-line trajectories as the 
molecules in the low field seeking state when no electric 



fields are applied, but are not included in the trajectory 
simulations. The ring diameters are accurately repro- 
duced by the simulations, but the simulated rings have 
sharper edges than are observed experimentally. Since 
the peaks at the ring radii are wider in the experimental 
data than in the simulated images, the contrast is also 
reduced. As a consequence some of the rings, in partic- 
ular the outermost ones, are barely visible in the exper- 
imental images while they are still clearly visible in the 
simulations. Simulations show that the rings stem from 
molecules that pass close to the electrodes at the end of 
the decelerator, and that are slightly focused. The in- 
tensity distribution in these rings depends sensitively on 
the electric field distribution just behind the decelerator. 
This region is more strongly influenced by the support 
structure of the electrodes, which is not included in the 
cylindrically-symmetric model. 

Overall, the measurements show only minor devia- 
tions from the predictions of the two-dimensional tra- 
jectory simulations, indicating that the misalignments of 
the decelerator are relatively small. Remaining differ- 
ences are attributed mostly to non-cylindrically symmet- 
ric fringe-fields behind the decelerator. These distort the 
transverse velocity distribution of the molecules emerg- 
ing from the decelerator by asymmetric beam-focusing or 
-defocusing. They have not been included in the present 
calculations, and to do so would indeed be quite demand- 
ing, both because the fields would have to be accurately 
known, and because the calculations would have to be 
done in three dimensions, without assuming cylindrical 
symmetry. 
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FIG. 6: Panel A: Angle-integrated radial intensity profiles 
of the two-dimensional transverse distribution at the detec- 
tor after actively guiding molecules at 288 m/s. The inset 
shows the narrow region below it with a higher time resolu- 
tion. Panel B: Results from trajectory simulations under the 
same conditions. 



be guided to the end of the decelerator. Panel B shows 
the results of numerical trajectory simulations that agree 
well with the experimental data. 



C. Guiding of molecules with AC fields 



D. Non-adiabatic losses 



Similar images to those shown in the previous section 
were also recorded after actively guiding the molecules at 
288 m/s. Panel A of figure [5] shows the angle-integrated 
radial intensity prohles of the transverse distribution for 
the range of arrival times covering the fast part of an ini- 
tial distribution centered around saSOO m/s. The main 
panel shows the distributions measured in 40 ps-wide ar- 
rival time gates over a large range. Such slices are ob- 
tained by applying a short voltage-pulse to the front of 
the MCP-stack in the imaging detector. The fast part 
of the initial molecular beam, which arrives around 2.2- 
2.3 ms, shows a narrow transverse distribution; these 
molecules move faster than the traps and are not confined 
by the longitudinal potential. They traverse the saddle 
points between traps and thus see a much weaker trans- 
verse confining potential. The inset of panel A shows the 
narrow range of arrival times below it with a higher res- 
olution (3 ps wide gates). The narrow peaks that result 
as the molecules are bunched while traversing the sad- 
dle points are visible here. The longitudinally-confined 
guided molecules are visible at an arrival time of 2.4 ms. 
Because these molecules are confined in a deep potential 
well, a much wider transverse velocity distribution can 



In the Stark decelerator on a chip, it was observed that 
non-adiabatic transitions from the low field seeking states 
to states only weakly affected by the electric fields result 
in significant losses of molecules from the traps. Similar 
losses have been observed for ammonia molecules trapped 
in macroscopic electrostatic traps, albeit on a longer time 
scale. "^^ The chip decelerator operates on a similar princi- 
ple as the one described here, in that the required fields 
for the deceleration are obtained via the application of 
sine-modulated potentials to a periodic electrode-array. 
It has been shown that, in order to decelerate molecules 
on this chip, the non-adiabatic losses must be prevented 
by either using the ^"^CO isotopologue^^ or by apply- 
ing an external magnetic field. "^^ Non-adiabatic losses of 
molecules from the traps on the chip result primarily from 
the small size of the traps: as the molecules pass near the 
zero-field trap center, the direction of the electric field 
vector changes rapidly, providing spectral components 
that can drive the transitions to non-trapped states. In 
the macroscopic traps of the decelerator presented here 
the rate of such transitions should be much lower. 

To confirm that the effect of such losses is small, 
i^C^^O and i^ci^Q, were both guided at 300 m/s and 
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FIG. 7: Comparison of the arrival time distributions for ^^CO 
and ^^CO for guiding at 300 ni/s (panel A) and deceleration 
from 300 m/s to 180 m/s in 2 ms (panel B). 



decelerated from 300 m/s to 180 m/s over the length 
of the decelerator, and the arrival time distribution of 
the molecules exiting the decelerator was recorded. Be- 
cause of the hyperfine structure in -'^^CO, any losses due 
to non-adiabatic transitions should be reduced in com- 
parison with 12C0,?6 

Figure [7lA shows, from bottom to top, the experimen- 
tal arrival time distribution for ^^CO and ^^CO, and the 
simulation for ^^CO when guiding of the molecules at 
300 m/s. In all three traces, the main peak shows the 
molecules that are moving in the main trap, at 300 m/s, 
and the sharp and regular structures to either side of the 
main peak are from molecules that are bunched as they 
cross the saddle points between the traps. Except for 
an overall intensity difference, the three traces show no 
significant differences. The source of the different intensi- 
ties lies mainly in different source conditions. Additional 
information on potential non-adiabatic losses can be ob- 
tained by comparing the traces for deceleration: this is 
shown in [TJB for ^'^CO (lower trace) and ^^CO (upper 
trace) . By normalizing these traces based on the relative 
intensities of the guiding data, it is found that about 60% 
more ^^CO is decelerated than ^'^CO. This difference in 
deceleration efficiency is qualitatively explained by the 
slightly lower mass of ^^CO. Significant non-adiabatic 
losses would lead to an enhanced efficiency of ^^CO de- 
celeration. 

Monte-Carlo trajectory simulations for guiding of 
^^CO (shown as uppermost trace of figure [TjA.) and ^^CO 
(not shown) show excellent agreement with the experi- 
mental data. The relative amounts of ^^CO and ^"^CO 
that can be decelerated also agree well with the pre- 
dictions of the simulations, which do not allow for non- 



adiabatic processes. On the timescale of the current ex- 
periments (which is similar to the radiative lifetime of 
the metastable molecular state that is used), losses due 
to non-adiabatic transitions are thus concluded to be neg- 
ligible for the case of CO. 



IV. SUMMARY AND OUTLOOK 

A detailed account on the design, construction, and op- 
eration of the recently developed traveling wave decelera- 
tor for polar neutral molecules has been given. The very 
high accuracy required for the mechanical design calls 
for special care during the production and alignment of 
the individual electrodes. Generating high-voltage sinu- 
soidal waveforms with frequencies between 25 kHz (for 
the current geometry) and ideally Hz is particularly 
challenging. In the currently employed design, the range 
12-26 kHz is accessible, and this could be extended by us- 
ing different transformers. To decelerate molecules from 
higher initial velocities or to lower final velocities one 
could segment the decelerator into different parts that are 
each driven by a separate set of transformers. By hav- 
ing slightly overlapping frequency ranges between neigh- 
boring segments, a smooth deceleration would then be 
possible over wider velocity ranges, albeit at the cost of 
flexibility. A more flexible solution, that however requires 
more development, will be the use of a linear amplifier 
that covers the entire range of interest. 

The good agreement between the simulations and ex- 
periments raises the confidence in the understanding of 
the overall dynamics of molecules inside this type of de- 
celerator. The longitudinal distribution is measured by 
recording the molecule arrival times, and the transverse 
distribution has been characterized via imaging of the 
complete packet of guided molecules in dc-electric fields. 
Narrow time-slices of the transverse distributions have 
been recorded for guiding in modulated fields, and it was 
found that the calculated and measured data agree rea- 
sonably well. Remaining discrepancies, for dc- and ac- 
guided molecules, are attributed to fringe fields that are 
generated from the electrode mounting structure. These 
fields are shielded from molecules that are inside the de- 
celerator, but not from molecules flying from the end 
of the decelerator to the detector. The comparison be- 
tween -'^^CO and ^■^CO showed that in the present geome- 
try, unlike the case of the chip decelerator, non-adiabatic 
losses can be completely neglected. This is because the 
increased trap-size leads to lower additional frequencies 
perceived by the molecules as they fly through the central 
region of the traps. 

As has been suggested previously,'^^ this type deceler- 
ator seems particularly suitable for the deceleration of 
heavy diatomic molecules. Each rotational state of these 
molecules has low- field- seeking components only below a 
certain electric field. Above that, all levels undergo an 
avoided crossing with levels from the next rotational state 
and become high-field-seekers, making them inaccessible 
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by traditional Stark deceleration. The only option so far 
was the use of alternating-gradient deceleration, but this 
approach has not yet been demonstrated to work stably 
and efhciently. In the present decelerator, however, the 
electric field magnitude inside the moving traps is always 
relatively low, yet the high electric field gradient never- 
theless enables efficient deceleration. 
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